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5 1 . 

o r Circumstellar disks are exposed to intense ultraviolet (UV) radiation from the young star. In 
q j the inner disks, the UV radiation can be enhanced by more than seven orders of magnitude 
fcj | compared with the average interstellar radiation field, resulting in a physical and chemical 
^ ' structure that resembles that of a dense photon-dominated region (PDR). This intense UV field 
affects the chemistry, the vertical structure of the disk, and the gas temperature, especially in 
the surface layers. The parameters which make disks different from more traditional PDRs are 
^ ! discussed, including the shape of the UV radiation field, grain growth, the absence of PAHs, 
OO \ the gas/dust ratio and the presence of inner holes. Illustrative infrared spectra from the Spitzer 
O Space Telescope are shown. New photodissociation cross sections for selected species, including 
simple ions, are presented. Also, a summary of cross sections at the Lyman a 1216 A line, 
known to be strong for some T Tauri stars, is made. Photodissociation and ionization rates are 



OG , computed for different radiation fields with color temperatures ranging from 30000 to 4000 K 
^ \ and grain sizes up to a few /im. The importance of a proper treatment of the photoprocesses 
^ ' is illustrated for the transitional disk toward HD 141569A which includes grain growth. 

X' 

- ■ 1 Introduction 

It is well established observationally and theoretically that many young stars in the solar 
neighborhood are surrounded by disks, i.e., flattened rotating concentrations of gas and dust 
(see [TS] for a review). In the early phases, up to a few Myr after collapse of the parent cloud, 



1 



the disks are rich in gas and dust inherited from the cloud core. Once accretion onto the 
star stops and the disk becomes less turbulent, the ~0.1 /im interstellar grains coagulate to 
larger and larger particles and settle to the midplane, eventually forming km-sized planetesimals 
which interact gravitationally to form protoplanets (e.g., [80]). In this period, between a few 
and 10 Myr, also most of the gas may be dissipated from the disk although the precise time 
scale for gas removal is still uncertain. Dust disks are observed around much older (up to a 
few Gyr) stars as well, but the grains in these disks are the 'debris' produced by collisions and 
fragmentation of planetesimals and are no longer the original interstellar particles (e.g., [37]). 
This paper is concerned with young, gas-rich disks up to the transitional stage when the disk 
becomes optically thin to UV radiation. 

Young stars are known to be powerful emitters of UV radiation, with typical intensities at the 
disk surface that are orders of magnitude higher than the average interstellar radiation field 
[23J. For some sources, most of the flux may be contained in the Lyman a line [3]. These UV 
photons play a very important role in the physical and chemical structure of the disk, especially 
in the surface and intermediate layers. Many disks are thought to have a flaring structure, in 
which the disk height H and opening angle H/R increase with radius R (e.g., [7J[TU]). In such a 
geometry the disk surface intercepts significantly more UV radiation from the star than in a flat 
geometry, heating the surface layers and increasing the scale height further when hydrostatic 
equilibrium is assumed in the vertical direction. Because the density in the upper layers is 
below that at which gas-grain coupling becomes effective, the gas temperature can be higher 
than the dust temperature due to photoelectric heating, affecting both the structure, chemistry 
and line formation J27J [2U]. In the inner disk (< 10AU), the temperature can even become so 
high (several thousand K) that photoevaporation becomes effective [17] . 

The main effect of the UV radiation on the chemistry is through photodissociation and pho- 
toionization of molecules and atoms in the surface layers. These processes remain effective until 
most photons have been absorbed by the dust and the photorates have fallen below those of 
other chemical reactions, typically at a depth of a few magnitudes of extinction. The chemical 
structure of the outer disk has been modeled by various groups in both the radial and vertical 
direction (e.g., [82l [T| WI\ [78] ) . In the vertical direction, the chemical structure resembles that of 
a photon-dominated region (PDR), with a transition from C + to C and CO around Ay wl mag 
and radicals such as CN and C2H having a significant abundance in the intermediate layers. In 
the cold mid-plane (< 20 K), most molecules are frozen out onto the grains. This layered chem- 
ical structure is consistent with observations of ions and radicals with abundance ratios that 
are similar or higher than those found in PDRs (e.g., [HJ [70]). If the UV radiation is energetic 
enough to dissociate CO, isotope selective processes may occur leading to fractionation of CO 
isotopes which can eventually be incorporated into meteorites |46j. Finally, UV photons can 
affect the chemistry through photodesorption of ices brought to the surface by vertical mixing 
in the disk [82]. 

The above discussion illustrates the need for accurate photodissociation cross sections to de- 
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scribe the disk chemistry. In this paper, photodissociation rates relevant for disks exposed 
to different types of radiation fields are presented, including calculations of cross sections for 
species not considered before. Also, a critical evaluation of available cross sections at the Ly- 
man a wavelength is given. The resulting chemical structure is illustrated for the transitional 
disk around the Herbig Ae/Be star HD 141569A. 

2 Photoprocesses in disks 

The PDR structure of a disk differs from that of the more commonly studied PDRs in molecular 
clouds [21] . First, the spectral shape of the stellar radiation field can differ significantly from 
that of the standard interstellar radiation field (ISRF), especially for the cooler T Tauri stars. 
Figure 1 compares the standard ISRF cf. [12] with that of a 10000 K (typical Herbig Ae star) 
and a 4000 K (typical T Tauri star) radiation field, normalized to have the same integrated 
intensity from 912-2050 A. Below 2000 A where most molecules are photodissociated, the 4000 
K field has orders of magnitude lower intensity. The figure also includes a NEXTGEN simulated 
spectrum of a B9.5 star [22], illustrating that the intensity may be further reduced by stellar 
absorption lines in the critical 912-1100 A where H 2 and CO are photodissociated. These 
stellar lines become increasingly important at lower stellar temperatures. Some T Tauri stars 
are known to have excess UV emission above that of a 4000 K blackbody, either from a hot 
boundary layer between the accretion disk and the star or from chromospheric stellar activity, 
bringing the overall shape back closer to that of the ISRF [6[ [30] . Specific resonance lines like 
Lyman a can dominate the spectrum [3]. The key chemical effect is that H 2 and CO are hardly 
photodissociated by 'cool' stars at 912-1100 A; only the general ISRF incident on the disk 
sets up the H— ► H 2 and C + ^C— >CO transitions [78J. Other molecules with photodissociation 
channels primarily at short wavelengths (e.g., N 2 , CN) are equally affected. 
Second, the UV field can be orders of magnitude higher than that studied in molecular clouds. 
For example, at 15 AU from an AO star, the radiation is 10 7 times that of the ISRF. For such 
high fluxes, photorates can become comparable to dissociative recombination rates so that 
photodissociation of ions needs to be included in the models. 

Third, the grain sizes and properties may be very different from those found in clouds. In 
particular, there is strong evidence from infrared observations that grains in the upper disk 
layers have grown from ~ 0.1 /mi to a few /im in size regardless of the spectral type of the star 
(e.g., [7H E2])- Figure 2 illustrates this growth using recent Spitzer Space Telescope spectra 
of the 10 and 20 /im silicate feature. Larger grains have significantly lower extinction at UV 
wavelengths. Also, PAHs may have lower abundances and are detected toward only few T Tauri 
stars [To] . 

Fourth, the gas/dust mass ratio can differ from the interstellar medium value of 100. Currently, 
it is unclear whether the (small) dust particles disappear before the gas or vice versa. However, 
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dust-free gaps, rings and holes in disks have been found for some objects including the HD 
141569A disk (e.g., [21 [8]). Some remaining gas can be present in these dust holes (e.g., |50j). 
Fifth, the density in the disk PDR varies by orders of magnitude, both in vertical and radial 
direction, in contrast with the near constant density assumed for clouds. If ices are released into 
the gas due to vertical mixing followed by thermal or photodesorption, the initial conditions for 
the PDR are different from those of standard clouds. For example, the presence of enhanced 
gaseous H2O has been shown to modify the PDR chemistry for the case of protostellar envelopes 

m 

Finally, young stars are known to be powerful X-ray emitters (e.g., [15J). In addition to direct 
ionization, the UV radiation resulting from the interaction of the secondary electrons with H2 
has a major effect on the chemistry [T9(, [63]. This spectrum consists of many discrete lines 
originating from the H2 B and C electronic states, with peaks at 1500-1600 A and shorter 
wavelengths. Thus, many of the basic photodissociation processes discussed here are also 
relevant for models which take X-rays into account. 



3 Photodissociation cross sections and rates 

The photodissociation rate of a molecule can be computed from 

*T = J a(X)I(X)d\ s^ 1 (1) 

where a is the photodissociation cross section in cm 2 and / is the mean intensity of the radiation 
in photons cm -2 s~ x A -1 . Usually a is a broad, continuous function with wavelength peaking 
close to the vertical excitation energy of the electronic state involved in the process. For 
photodissociation initiated by line absorptions (e.g., predissociation), the rate becomes 

71 f 2 

k l ;r = —>>lifuiVuiut s- 1 (2) 

' mc z 

where f u £ is the oscillator strength for absorption from lower level t to upper level u and r\ u the 
dissociation efficiency of state u, which lies between and 1. The numerical value of the factor 
Tie 2 /mc 2 is 8.85 x 10~ 21 in the adopted units with A in A. The total photodissociation rate of 
a molecule is obtained by summing over all channels. 

Overviews of photodissociation cross sections and interstellar photodissociation rates of astro- 
physically relevant molecules have been given by [101 [73], E7J [26] . Data on cross sections can be 
found in the chemical physics literature, either from experiments (stable molecules) or theory 
(radicals, ions). These summaries include only limited data on the photodissociation of small 
ions which are key to building carbon- and oxygen species. Therefore, a literature search was 
performed and new cross sections were calculated for several species. Also, the first critical 
evaluation of cross sections at the Lyman a wavelength is presented. Finally, photodissociation 
rates are calculated for different radiation fields and grain parameters. 
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3.1 New photodissociation cross sections 

For small molecules and ions, quantum chemical calculations of the potential energy curves 
and transition dipole moments combined with dynamical calculations of the nuclear motions 
can provide accurate photodissociation cross sections and oscillator strengths (see [M] for re- 
view). All new calculations presented here were performed with the MOLPRO set of programs 
[81], using the VTZ atomic orbital basis set [13] . For neutral molecules, diffuse s and p func- 
tions were added to allow for a proper description of molecular Rydberg states. Molecular 
orbitals were generated in state-averaged complete active space (CAS) calculations. In the 
ultimate contracted multi-reference configuration-interaction (CI) calculation, typically 50 ref- 
erence configurations per symmetry were used. Thus, around 100000 contracted configurations 
out of 300000 uncontracted configurations were generated per symmetry. Orbitals with an 
orbital energy below -15 eV were kept doubly occupied in all CI calculations. 
Typically, the lowest 5 electronic states of each symmetry were calculated at the ground-state 
equilibrium geometry together with the corresponding transition moments. For diatomics, 
complete potential energy curves were computed as well. Many of the excited states are bound 
and the predisssociation efficiencies are often unknown. In the following, individual cases are 
discussed in more detail. 

CH^: Accurate CI calculations by Theodorakopoulos & Petsalakis [69] show many dipole- 
allowed excited states below 13.6 eV. Their 1, 2 and 3 2 B 2 , 2, 3 and 4 2 Ai and 2 2 Bx states have 
been included with /=0.008, 0.0001, 0.02, 0.01, 0.06, 0.05 and 0.03, respectively. 
CH4 : Detailed studies of the CH4 photodissociation processes starting from its lowest C 2v 
state have been carried out by van Dishoeck et al. [76] indicating that the higher excited states 
are likely dissociative. The oscillator strengths into the 2, 3 2 Ai and 2 2 Bi states computed in 
this work are /=0.04, 0.04 and 0.08, respectively. 

O^: Calculations show that all electronic states below 13.6 eV are bound [2S|. Our computed 
oscillator strengths to the A 2 II n and 2 2 II U states — the only dipole allowed states below 13.6 
eV and above the dissociation limit- are only /=0.005 each. 

CO + : CO + is strongly bound with a dissociation energy of 8.34 eV, but its higher excited D 2 I1, 
G, E and F 2 £ + states are likely (pre-)dissociative [38]. The calculated oscillator strengths to 
the D and G states are /=0.01 and 0.02, respectively. Similar values have been assumed for 
the higher states. 

HCO + : The excited electronic states and photodissociation processes of HCO + have been 
studied in detail by Koch et al. [351 EH]. HCO + is remarkably transparent at UV wavelengths: 
the only dipole- allowed dissociative state is the l 1 !! state around 11.5 eV with a small cross 
section. 

H 2 + : Calculations show that there are no dipole-allowed dissociative electronic states below 

13.6 eV so that the interstellar H 2 + photodissociation rate is negligible. 

SiO: The electronic structure of SiO is similar to that of CO, but with a lower dissociation 
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energy of 8.26 eV implying that even the lower Rydberg states can contribute to SiO photodis- 
sociation if they are fully predissociated. Our computed oscillator strengths to the 3 1 X + and 
2, 3, 4 and 5 x n states are /=0.10, 0.32, 0.03, 0.11 and 0.10, respectively. 

3.2 Lyman a cross sections 

The Lyman a line can outshine the continuum radiation, so that it is important to know 
whether a molecule can be photodissociated at 1216 A or not. Table 1 summarizes the cross 
sections for a number of key molecules, including those that are known to have significant 
abundances in the inner disk [47J. The table contains references to the original experimental 
or theoretical work as well as an assessment of the accuracy of the data, with the notation A 
(cr known to better than 50%), B (cr uncertain up to a factor of two) or C (a uncertain up to 
an order of magnitude). 

The following molecules cannot be dissociated by Lyman a radiation: H2, CO, N2, CN, H^, 
OH + , H 2 + and HCO + . An interesting case is formed by O2. The absorption cross section at 
1215-1216 A is known to be small, of order (1 — 2) x 10~ 20 cm 2 [53], but there are resonances 
into the v' = and 1 levels of the 2 3 £~ state at 1245 and 1205 A, with / = 0.015 each [5]. 
If the Lyman a line is as wide as observed for some stars (1205-1230 A, ±0.1 eV [3]), these 
states can make a small contribution. 

Experimentally determined cross sections at Lyman a are generally accurate to better than 
a factor of 2, especially when the absorption is continuous. For molecules for which only 
theoretical calculations are available (e.g., CH, CH 2 , C2H, C3, ...), the cross sections at Lyman 
a are highly uncertain because the energy level calculations are uncertain by ~0.2 eV, which is 
not accurate enough to determine whether there is an exact resonance of a predissociative state 
at 1216 A. For continuous absorption into a directly repulsive state around 10.2 eV (e.g., OH), 
the computed results are reliable, however. For some molecules, no estimate is possible because 
of lack of both experimental and theoretical data around Lyman a. This includes C3H2, SH, 
SH+, CS and SiO. 

Only a few molecules can be photoionized by Lyman a (see Table 2 of [73] for ionization 
potentials). For species like NH3 and NO, the photoionization cross sections or yields have 
been measured so that the dissociation and ionization parts can be separated. In other cases 
such as CH3, the branching ratio is unknown. SO has an ionization potential of 10.29 eV, just 
above the Lyman a threshold, whereas CS2 has 10.08 eV, just below. Of the major atoms, C, 
N, O, and S cannot be photoionized by Lyman a whereas Mg, Si and Fe can. 

3.3 Photodissociation rates for Herbig Ae and T Tauri stars 

The database on photodissociation and ionization cross sections assembled by van Dishoeck 
[73] has been used together with the above new input to compute photodissociation rates for 
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various radiation fields. No attempt has been made to systematically review the chemical 
physics literature since 1988 for molecules other than those discussed in §3.1. Minor updates 
are included for some species such as CH 2 , for which the cross sections of [77J are adopted. 
References to the sources for the cross sections are given in Table 1 of [73] , which also includes 
the most likely dissociation products. The resulting rates for the standard ISRF cf. 1 2] are 
given in Tables 2 and 3. 

In calculating the photodissociation rates, it has been assumed that all states above the dis- 
sociation limit have 100% dissociation efficiency, thus providing a maximum value (rj u = 1). 
Even with this assumption, the photodissociation of ions is on average significantly slower than 
that of neutrals. Several ions (e.g., FLt, H 2 + ) do not even have any dipole-allowed dissocia- 
tive electronic transitions below 13.6 eV. Note that CH + dissociates primarily into C + H + 
for the radiation fields considered here, not C + + H as listed in UMIST99. In contrast, OH + 
dissociates primarily into + + H. 

The photodissociation rates for larger molecules are uncertain because experimental data are 
sparse and incomplete. Moreover, often only absorption cross sections are known, not the 
dissociation yields. For large molecules, absorption is followed by internal conversion to highly 
excited vibrational levels of the ground electronic state, only some of which lead to dissociation. 
Also, the products of the photodissociation are largely unknown and can vary with wavelength; 
see [48] for an illustrative example for the case of CH3OH. 

Tables 2 and 3 include the photodissociation and photoionization rates for a Tbb — 10000 K 
and a 4000 K blackbody radiation field. As noted in §2, the shapes of the UV fields of Herbig 
Ae and T Tauri stars are thought to lie in between these three extremes. Their intensities have 
been normalized such that the integrated values from 912-2050 A are the same as those of the 
Draine [T2j field, 2.67 x 10~ 3 erg cm -2 s _1 . The latter value is a factor of 1.7 larger than that 
of the integrated Habing [21] field of 1.6 x 10~ 3 erg cm -2 s _1 used in other normalizations. The 
adopted dilution factors are 1.57 x 10~ 14 and 1.90 x 10~ 9 for 10000 and 4000 K, respectively. 
For applications in disk chemistry, these rates need to be scaled to the appropriate strength at 
a certain distance from the star. The photorates for (scaled) fields with Tbb = 20000 — 30000 K 
are close to those for the ISRF. 

As shown previously for interstellar clouds exposed to cool stars [61] , the photorates can decrease 
by orders of magnitude if a 4000 K blackbody is used. Species like H 2 , CO, C, CN and N 2 , which 
are only dissociated or ionized at <1100 A, are most affected (cf. Figure 1): their photorates 
become negligible for Tbb = 4000 K. In contrast, molecules with dissociation channels at 
wavelengths as long as 3000 A (e.g., CH, HCO, O3) can have enhanced rates in the cooler fields 
with the adopted normalization. Molecules which absorb over a broad wavelength range (e.g., 
OH, H 2 0, NO) vary by only a factor of a few between Tbb — 30000 and 4000 K. 
Illustrative examples of the effect of Lyman a radiation on the photorates in disks are given by 
[3J. The photodissociation and ionization efficiencies corresponding to Lyman a absorption in 
the context of X-ray induced chemistry are summarized by [H] for selected species. 
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3.4 Depth-dependence of photorates 

With depth into a cloud or disk, the UV radiation is attenuated due to absorption and scattering 
by dust grains [56]. For typical 0.1 /im interstellar grains, the extinction, albedo and scattering 
phase function are taken from [57J . The rates are then fitted to a single exponential decay with 

Kd = ^exp(-7Ay) s _1 . (3) 

The fits of 7 are performed over the range Ay = — 3 mag and are appropriate for a model 
in which the radiation is incident from one side only. The absolute values of 7 vary by ~ 10% 
depending on the extinction range chosen for the fits; the relative values from molecule to 
molecule are more accurate. Tables 2 and 3 contain the results for the different radiation fields. 
The values of 7 for the ISRF differ from those of van Dishoeck [73] due to the adopted grain 
properties, which are in between those of grain models 2 and 3 of [56] used previously. Inter- 
estingly, the values of 7 do not change dramatically with T BB for most species, even though 
there is a general trend for 7 to become smaller with cooler fields (see also [61J). 
When grains grow larger, the main effect on the grain parameters is the smaller extinction at 
UV wavelengths, which becomes comparable to that at visible wavelengths (e.g., [59]). The 
photorates have been re-computed as a function of depth using the properties for /zm-sized 
ice-coated grains derived for the HD 141569A disk [13, [2SJ. For species such as C and CO 
which absorb only at the shortest wavelengths, the exponents 7 are lowered from > 3 to ~0.6 
for the ISRF. For molecules which are dissociated primarily at long wavelengths (e.g., HCO), 
7 is lowered from 1.1 to 0.47. Other species have exponents 7 in between these two values. 
For Tbb — 4000 K, the trends are similar, with values of 7 between 0.6 and 0.38. Thus, the 
differences between molecules in the depth dependence of the photorates are minimized for 
large grains. 

4 A disk model example 

In this work, the effects of different radiation fields and different grain properties are taken into 
account explicitly by calculating at each point in the disk the UV radiation field from the star 
as well as the ISRF (either in 1+1D or in full 2D) and then computing the photodissociation 
rates by integrating the cross sections multiplied by the radiation intensity over wavelength 
[751 125] . Self-shielding of H 2 and CO as well as mutual shielding of CO and its isotopes are 
included. X-rays are not included, nor is time- dependence, vertical mixing or photodesorption 
of ices. 

As an example, the disk around the Herbig Ae/Be star HD 141569A (B9.5 V, 22 L , 99 pc) 
is modelled. This ~5 Myr old star has a large ~500 AU disk with a huge inner hole in the 
dust out to 150 AU and two dust rings at 185 and 325 AU seen in scattered light images (e.g., 
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[2, 8_p). The outer rings may be explained by tidal interactions with two M-type companions, 
HD 141569B and C, but other explanations such as a giant planet are not excluded. The origin 
of the inner hole is unknown. The disk is optically thin to UV continuum radiation and the 
grains have grown to at least /im size [45J . Nevertheless, gas is still present since the CO J=2- 
1 and 3-2 millimeter [TT] and the v — 1 — infrared [3] lines have been detected. Moreover, 
PAHs have been seen, albeit at a low level [SB]- Thus, this disk provides a good opportunity 
to investigate the effects of different UV radiation fields and larger grains. 
Detailed models to constrain the total gas mass in this disk have been developed by Jonkheid et 
al. [28]. They use a 1+1D approach, with ID PDRs computed in the radial rather than vertical 
direction for this optically thin disk. Both the radiation from the star (see Figure 1) and the 
ISRF are included. The size of the dust grains enters the chemistry in various ways: in the UV 
extinction, in the H 2 formation rate, and through the photoelectric heating efficiency [31]. The 
best-fitting model to the CO millimeter lines has a total gas mass of 80 Msarth, compared with 
a total dust mass of 2.2 MEarth in grains up to 1 cm. An important conclusion is that some gas 
must be present in the inner hole to provide sufficient H 2 and CO self-shielding in the radial 
direction. The PAH abundance with respect to total hydrogen is ~ 1CT 10 , but even at this low 
abundance PAHs are an important site for H 2 formation. Models without PAHs require much 
larger gas masses [28] . 

Figures 3 and 4 present the radial and vertical distributions of various molecules in the HD 
141569A disk, both with and without the ISRF. Because the stellar radiation field has few 
photons at wavelengths below 1200 A, the photodissociation of CO and H 2 and photoionization 
of C occurs mostly by the ISRF in the outer disk. The vertical slices show the typical PDR 
structure with atomic H and C + dominant on the surface, and the transition to H 2 taking place 
in the intermediate layer. However, the CO column never becomes large enough in the vertical 
direction for significant self- shielding so that carbon stays in atomic form. Other molecules 
like CH and C 2 H which are not self-shielding and can photodissociate over a large wavelength 
range continue to be dissociated very rapidly by the stellar radiation. Thus, even without the 
ISRF, the CO abundance stays low since its precursor molecules have very low abundances. 
Neutral atomic carbon is the dominant form of carbon since there are few ionizing photons. The 
predicted peak line intensity of ~1 K should be readily detectable by submillimeter telescopes 
[28J. Failure to detect this line may imply that the observed CO is not a relic of the interstellar 
cloud from which HD 141569A formed. 



5 Conclusions 

The main results of this work are as follows. 

*see http://hubblesite.org/newscenter/newsdesk/archive/releases/2003/02/ for image 
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• New photodissociation cross sections are presented for several species not considered 
before. Also, the first critical evaluation of cross sections at Lyman a is given. 

• The photodissociation of ions is potentially important in (inner) disk chemistry where 
the radiation field can be 10 7 times more intense than the standard ISRF. However, 
the photorates of many ions, especially those containing oxygen, are found to be slow 
compared with those of neutrals. 

• Photorates have been computed for a range of radiation fields appropriate for Herbig 
Ae/Be and T Tauri stars. Molecules such as H 2 , CO, N 2 and CN, which are photodisso- 
ciated only below 1200 A, have rates that are more than 5 orders of magnitude decreased 
for a T BB = 4000 K blackbody field. 

• The exponent 7 characterizing the depth-dependence of the photorates decreases with 
Tbb but the effect is not large in the 4000-30000 K range for most species. For larger 
//m-sized grains, the depth dependence of the photorates becomes much shallower, with 
7 falling in a narrow range of 0.4-0.6 for all species. 

• The example of the HD 141569A disk illustrates that both the shape of the stellar radia- 
tion field and the size of the grains affect the chemistry. Both need to be treated correctly 
to derive quantitative conclusions about the gas mass and chemistry from observed lines. 
Even a small amount of PAHs or small grains can significantly affect dissociation and 
ionization rates, as well as the H 2 formation rate. 
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Table 1: Photodissociation and ionization cross sections at Lyman a 1216 



Species 


(cm -2 ) 


Accuracy 6 


Ref 




CH 


5.0(-20) 


C 


1721 


T 


CH 2 


5.0(-20) 


C 


[771 


T 


CH 4 


1.8(-17) 


A 


[4T1 
1 i 


E 


c 2 


5.0(-18) 


B 


[55] 


T 


c 3 


1.0(-18) 


C 


[58] 


T 


C 2 H 


1.0(-18) 


C 


1601 

I i 


T 


C2H2 


>4(-17) c 


B 


1651 


E 


C4H2 


3.5(-17) 


B 


1541 


E 


OH 


1.8(-18) 


B 


[75\ 


T 


H 2 


1.2(-17) 


A 


1421 

L J 


E 


2 


1.0(-20) b 


C 


[53J 


E 


C0 2 


6.1(-20) 


A 


[83J 


E 


H 2 CO 


1.0(-17) 


B 


1671 


E 


CH3OH 


1.4(-17) 


A 


[52] 


E 


NH 


1.0(-18) 


B 


[33] 


T 


NH 3 


1.0(-17) 


A 


[M] 


E 


HCN 


3.0(-17) 


A 


[39] 


E 


HC 3 N 


2.5(-17) 


B 




E 


CH 3 CN 


2.0(-17) 


A 


[66J 


E 


NO 


4.0(-19) 


B 


m 


E 


H 2 S 


3.3(-17) 


B 


[40\ 


E 


SO 


1.0(-16) 


C 


[51] 


E 


S0 2 


3.0(-17) 


B 


[40J 


E 


ocs 


1.5(-17) 


A 


[40J 


E 


cs 2 


2.5(-17) 


B 


1 A r\l 

m 


E 


Mg p.i. 


3.0(-19) 


A 


[75] 




Si p.i. 


3.0(-17) 


A 


[73] 




Fe p.i. 


6.2(-19) 


A 


[73] 




NH 3 p.i. 


2.0(-18) 


B 


[79] 


E 


NO p.i. 


1.6(-18) 


B 


[20] 


E 


CS 2 p.i. 


2.0(-16) 


B 


[40j 


E 



See text for list of species that cannot be dissociated or ionized by Lyman a 

See text for discussion 

/=0.013 
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Table 2: Photodissociation rates for various radiation fields"' 1 



Species 






K pd 


(s- 1 ) 








7 




ISRF C 


10000 K d 


4000 K d 


ISRF 


10000 K 


4000 K 


H+ 


5.7{ 


-10) 


1.9( 


-10) 


2.9( 


-11) 


2.37 


2.14 


1.99 


CH 


9.2( 


-10) 


2.0( 


-9) 


1.2( 


-7) 


1.72 


1.49 


1.28 


CH+ 


3.3( 


-10) 


3.5( 


-11) 


4.8( 


-10) 


2.94 


1.78 


1.31 


CH 2 


5.8( 


-10) 


1.2( 


-9) 


2.1( 


-9) 


2.02 


2.02 


2.12 


CH+ 


1.4{ 


-10) 


7.4( 


-11) 


2.6( 


-11) 


2.21 


1.91 


1.88 


CH 3 


2.7{ 


-10) 


2.5( 


-10) 


8.2( 


-10) 


2.27 


2.24 


2.32 


CH 4 


1.2{ 


-9) 


2.2( 


-10) 


1.2{ 


-12) 


2.59 


2.45 


2.29 


CH+ 


2.8( 


-10) 


4.2( 


-11) 


1.3( 


-13) 


2.71 


2.58 


2.48 


c 2 


2.4( 


-10) 


4.1( 


-11) 


3.2( 


-13) 


2.57 


2.36 


2.25 


C 2 H 


5.2( 


-10) 


1.9( 


-10) 


7.2{ 


-12) 


2.30 


2.16 


2.10 


C 2 H 2 


3.3( 


-9) 


1.2( 


-9) 


1.3{ 


-10) 


2.27 


2.12 


1.97 


C 2 H 4 


3.0( 


-9) 


2.2( 


-9) 


5.2( 


-10) 


2.10 


1.96 


1.90 


c 3 


3.8( 


-9) 


2.9( 


-9) 


2.0( 


-10) 


2.08 


2.07 


2.06 


c-C 3 H 2 


1.9{ 


-9) 


1.7( 


-9) 


9.2( 


-10) 


2.07 


2.06 


2.10 


OH 


3.9( 


-10) 


1.8( 


-10) 


1.3{ 


-10) 


2.24 


2.00 


1.67 


OH+ 


1.1( 


-11) 


7.8( 


-13) 


5.8( 


-13) 


3.50 


2.80 


1.75 


H 2 


8.0( 


-10) 


4.3( 


-10) 


1.2( 


-10) 


2.20 


1.97 


1.90 


o 2 


7.9{ 


-10) 


4.9( 


-10) 


4.5( 


-11) 


2.13 


2.05 


1.97 


o 2 + 


3.5( 


-11) 


3.6( 


-11) 


1.0( 


-11) 


2.02 


1.92 


1.90 


H0 2 


6.7{ 


-10) 


1.9( 


-9) 


1.2( 


-8) 


2.12 


2.08 


1.99 


H 2 2 


9.5( 


-10) 


4.3( 


-10) 


1.4( 


-10) 


2.28 


2.07 


1.97 


3 


1.9( 


-9) 


5.4( 


-9) 


1.5( 


-7) 


1.85 


1.69 


1.57 


CO 


2.0( 


-10) 


1.5( 


-11) 


1.4{ 


-15) 


3.53 


3.47 


3.24 


CO+ 


1.0( 


-10) 


2.2( 


-11) 


1.2{ 


-13) 


2.52 


2.43 


2.32 


co 2 


8.9( 


-10) 


9.0( 


-11) 


1.2( 


-12) 


3.00 


2.53 


2.00 


HCO 


1.1( 


-9) 


2.5( 


-9) 


3.5( 


-6) 


1.09 


1.09 


0.81 


HCO+ 


5.4( 


-12) 


4.5( 


-13) 


7.9( 


-17) 


3.32 


3.32 


3.32 


H 2 CO 


1.0{ 


-9) 


6.7( 


-10) 


1.8( 


-10) 


2.16 


1.99 


1.90 


CH3OH 


1.4{ 


-9) 


5.9( 


-10) 


7.0( 


-11) 


2.28 


2.07 


1.95 
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Table 3: Photodissociation rates for various radiation fields 0- ' 6 (Table 2 cont'd) 



Species 








(s- 1 ) 








7 






ISRF C 


10000 K d 


4000 K d 


ISRF 


10000 K 


4000 K 


NH 


5.0( 


-10) 


1.6( 


-10) 


3.0( 


-12) 


2.33 


2.24 


2.12 


NH+ 


5.4( 


-11) 


1.8( 


-10) 


8.4( 


-9) 


1.64 


1.52 


1.51 


NH 2 


7.5( 


-10) 


1.0( 


-9) 


5.7( 


-10) 


2.00 


1.90 


1.89 


NH 3 


1.2( 


-9) 


1.0( 


-9) 


1.1( 


-9) 


2.12 


1.99 


2.00 


N 2 


2.3( 


-10) 


1.4( 


-11) 


3.0( 


-16) 


3.88 


3.89 


3.87 


NO 


4.7( 


-10) 


4.3( 


-10) 


2.9( 


-10) 


2.12 


1.96 


1.94 


N0 2 


1.4( 


-9) 


1.0( 


-9) 


3.4( 


-10) 


2.12 


1.97 


1.92 


N 2 


1.9( 


-9) 


4.8( 


-10) 


2.0( 


-11) 


2.44 


2.32 


2.02 


CN 


2.9( 


-10) 


2.1( 


-11) 


2.0( 


-15) 


3.54 


3.49 


3.23 


HCN 


1.6( 


-9) 


2.5( 


-10) 


3.7( 


-12) 


2.69 


2.44 


2.02 


HC 3 N 


5.6( 


-9) 


3.0( 


-9) 


2.5( 


-10) 


2.16 


2.12 


2.12 


CH 3 CN 


2.5( 


-9) 


4.8( 


-10) 


8.5( 


-12) 


2.58 


2.38 


2.01 


SH 


9.8( 


-10) 


1.3( 


-9) 


1.6( 


-8) 


2.04 


1.85 


1.34 


SH+ 


2.5( 


-10) 


3.4( 


-10) 


4.0( 


-8) 


1.66 


1.30 


1.29 


H 2 S 


3.1( 


-9) 


2.0( 


-9) 


3.2( 


-9) 


2.27 


2.12 


2.16 


cs 


9.8( 


-10) 


2.7( 


-10) 


3.7( 


-12) 


2.43 


2.33 


2.14 


cs 2 


6.1( 


-9) 


1.3( 


-8) 


2.0( 


-8) 


2.06 


2.02 


2.03 


ocs 


3.7( 


-9) 


3.1( 


-9) 


7.0( 


-10) 


2.07 


1.98 


1.94 


so 


4.2( 


-9) 


4.4( 


-9) 


9.4( 


-9) 


2.37 


2.18 


2.16 


so 2 


1.9( 


-9) 


7.4( 


-10) 


2.8( 


-10) 


2.38 


2.11 


1.94 


SiH 


2.8( 


-9) 


1.4( 


-8) 


8.0( 


-7) 


1.59 


1.55 


1.22 


SiH+ 


2.7( 


-9) 


1.0( 


-8) 


3.3( 


-6) 


1.21 


1.13 


1.11 


SiO 


1.6( 


-9) 


5.6( 


-10) 


9.9( 


-12) 


2.28 


2.21 


2.19 



See van Dishoeck [73] for products and references to cross section data 

H^, HeH + and H 2 + cannot be photodissociated by radiation with A > 912 A, so k° t 

ISRF according to [12] with extension at >2000 A of [7J] 

Scaled blackbody radiation field with temperature Tbb (see text) 
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Table 4: Photoionization rates for various radiation fields" (Table 3) 



Species 








(O 








7 




ISRF" 


10000 K c 


4000 K c 


ISRF 


10000 K 


4000 K 


C 


3.1( 


-10) 


2.5( 


-11) 


4.2( 


-15) 


3.33 


3.27 


3.10 


Mg 


7.9( 


-11) 


5.9( 


-11) 


6.9( 


-12) 


2.08 


2.00 


1.96 


Si 


3.1( 


-9) 


1.2( 


-9) 


4.1( 


-11) 


2.27 


2.17 


2.09 


S 


6.0( 


-10) 


5.9( 


-11) 


4.0( 


-14) 


3.08 


2.95 


2.76 


Fc 


2.8( 


-10) 


1.3( 


-10) 


5.8( 


-12) 


2.20 


2.14 


2.05 


CH 


7 fil 


1(Ti 




1 1 "1 
-111 


l.U 




3 28 


3 20 


2 97 


CFL 


u.o 


19"l 
-1ZJ 




1 1\ 


o.u 




3 94 


3 94 


3 93 


Co 


A 1 1 


-1UJ 




-111 


l.U 


-10) 


3 81 


3 81 


3 78 


C2H2 


3.3( 


-10) 


2.4i 


-11) 


2.4( 


-15) 


3 52 


3 49 


3 33 


C2H4 


4.1( 


-10) 


3.5( 


-11) 


1.1( 


-14) 


3.21 


3 1 1 


2.91 


C2H6 


2.3( 


-10) 


1.5( 


-11) 


7.6( 


-16) 


3.74 


3.73 


3.55 


2 


7.6( 


-11) 


4.8( 


-12) 


1.2( 


-16) 


3.87 


3.87 


3.85 


H 2 


3.1( 


-11) 


2.0( 


-12) 


3.6( 


-17) 


3.90 


3.90 


3.88 


NH 3 


2.8( 


-10) 


2.6( 


-11) 


1.1( 


-14) 


3.12 


3.04 


2.86 


NO 


2.6( 


-10) 


2.9( 


-11) 


5.4( 


-14) 


2.93 


2.71 


2.42 


N0 2 


1.5( 


-10) 


1.2( 


-11) 


2.9( 


-15) 


3.33 


3.22 


2.87 


N 2 


1.7( 


-10) 


1.1( 


-11) 


1.4( 


-16) 


3.93 


3.93 


3.92 


H 2 S 


7.3( 


-10) 


7.2( 


-11) 


3.5( 


-14) 


3.09 


3.01 


2.86 


cs 2 


1.7( 


-9) 


1.5( 


-10) 


8.3( 


-14) 


3.16 


3.02 


2.77 


ocs 


6.9( 


-10) 


1.4( 


-11) 


8.5( 


-15) 


3.35 


3.29 


3.10 


H 2 co 


4.8( 


-10) 


4.1( 


-11) 


1.2( 


-14) 


3.21 


3.13 


2.96 



a See van Dishoeck |73j for references to cross section data 
b See footnote c Table 2 
c See footnote d Table 2 
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Fig. 1 Comparison of the interstellar radiation field (ISRF) according to [12] (with the 
extension by [74] for >2000 A) with a 10000 K (dashed) and 4000 K (dotted) black- 
body scaled to have the same integrated intensity from 912-2050 A. The scaled 
NEXTGEN model radiation field a B9.5 star [22J (as appropriate for HD 141569A) 
is included as well (dash-dotted). 

Fig. 2 Evidence for grain growth in disks around T Tauri stars. Top: Spitzer Space Tele- 
scope observations of the 10 and 20 fim silicate Si-0 stretching and O-Si-0 bending 
mode features of two T Tauri stars, shifted by +0.4 and +0.2 for clarity. Bottom: 
normalized absorption efficiencies Qahs for models of spherical amorphous olivines 
with various sizes, calculated using the distribution of hollow spheres (DHS) method 
of [49]. The models are offset by +0.4, +0.25 and +0.1, respectively. Figure based 
on [32]. 

Fig. 3 Left: Radial distribution in the midplane of the HD 141569A disk of the tempera- 
ture and density (top), and various chemical species (middle and bottom). Right: 
Vertical slice at R = 300 AU. This model includes both the stellar radiation (Figure 
1) and the ISRF. 

Fig. 4 As Figure 3, but without the ISRF 



19 



io L 



10' 



10 6 I 




h 2 ,co 

C p.i. 



CH 



HCO 



10000 K 



10" 5 ! ! 

i' i '/ 

■ \ ■ .•■ 

i ■ i .•• 

i 1 i / 

1(T lLl^l 



ISRF 


^ - ' ^-"-^ 












i ■ vim 
> 'I'l 




1 1/ 


! >; m ; 1 i 


'j ii! 


i i ' 

i ii / 




i ' ! 


i ii! 

ii 
1 1 



-4000 K 



1000 



2000 
Wavelength (A) 



3000 



Figure 1: Comparison of the interstellar radiation field (ISRF) according to [12] (with the 
extension by [74] for >2000 A) with a 10000 K (dashed) and 4000 K (dotted) blackbody scaled 
to have the same integrated intensity from 912-2050 A. The scaled NEXTGEN model radiation 
field a B9.5 star [22] (as appropriate for HD 141569A) is included as well (dash-dotted). 
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Figure 2: Evidence for grain growth in disks around T Tauri stars. Top: Spitzer Space Tele- 
scope observations of the 10 and 20 /im silicate Si-0 stretching and O-Si-0 bending mode 
features of two T Tauri stars, shifted by +0.4 and +0.2 for clarity. Bottom: normalized absorp- 
tion efficiencies Q a bs for models of spherical amorphous olivines with various sizes, calculated 
using the distribution of hollow spheres (DHS) method of [49] . The models are offset by +0.4, 
+0.25 and +0.1, respectively. Figure based on [32] . 
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Figure 3: Left: Radial distribution in the midplane of the HD 141569A disk of the temperature 
and density (top), and various chemical species (middle and bottom). Right: Vertical slice at 
R = 300 AU. This model includes both the stellar radiation (Figure 1) and the ISRF. 
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Figure 4: As Figure 3, but without the ISRF. 
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